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U  is  asymptotically  normal,  gamma  or  extreme  value  distributed.  The 
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decomposition  is  employed  to  provide  an  algorithm  for  generating  random 

8  distributed  numbers.  Many  of  the  results  are  based  on  insights 
/v/v ,  w 

provided  by  the  classical  theory  of  the  Gamma  function  in  the  complex 
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50.  Introduction  and  summary 

The  Beta  distribution  plays  a  key  role  in  multivariate  analysis  [2], 

fi27~ and  in  order  statistics- [14] .  A  useful  tool  for  the  asymptotic  study 

of  the  beta  variate  6  is  its  logarithm  U  *  -log  8  which,  as  we 

f  W  ~V  |  W  |  w 

will  see,  has  simple  structural  properties.  The  beta  variate  B  has 

,  w 

p.d.f . 

(0-1)  f£;v,wW  =  bTvTwT  xV~1U-x)W~1  ,  0  <  x  <  1  ,  0  <  v,w  , 


where  B(v,w)  is  the  beta  function  B(v,w)  =  r(v)r(w)/r(v+w) .  Correspondingly, 

the  density  for  U  is 
* » w 


(0.2) 


fU;v,«W 


e%v.w"'X> 


-vx 


u  -  > 


-XjW-1 


B(v,w) 


0  <  x  <  »  ,  0  <  v,w  . 


The  variate  JJ^  ^  has  the  simple  generating  function  $u>v  w(s) 

E [j6S]  *  /  fft(x)xsdx,  i.e., 

0  & 


-sU 

E[e  ~]  = 


(0.3) 


*JJ;v,w(s) 


r(v+w)  r(v+s) 
r(v)  r(v+w+s)  * 


Re  (s)  >  -v 


The  transform  (0.3)  is  the  basis  for  the  asymptotic  and  structural  study 
which  follows.  We  will  see  that  y  and  hence  £  have  a  simple  decomposition. 
One  finds  that 

(0.4)  y  «  yCM  ♦  yPF  , 

CM  PF  "CM 

whete  y  and  y  are  independent.  When  y  *  0,  it  has  a  completely 


»• 


v V 


■  « 


n 


t-' 
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PF 

monotone  p.d.f.  When  U  *  0,  it  has  a  p.d.f.  which  is  PF  in  the  nota- 
tion  of  total  positivity  [4].  On  the  basis  of  this  decomposition,  one 
sees  that  U  is  infinitely  divisible,  and  J3  infinitely  factorable  in  the 
corresponding  sense. 

The  powerful  apparatus  of  the  Gamma  function  in  the  complex  plane 
permits  one  to  find  the  asymptotic  behavior  of  ^  and  ^  as  v  and 
w  go  to  infinity.  The  behavior  is  simple  and  interesting.  It  will 
be  shown  that 

d  „  ...  ,  - 


(a) 

U 

~vw 

-  log 

(b) 

vU 

~v,w0 

*  *w 

1 

(c) 

2  - 

U 

~v,w 

~v,w 

0, 

N(0,1)  as  v,w  -*•  for  a  broad  simple 


U;v,w 

family  of  paths  given  in  Section  2. 

In  (a),  G  is  a  conjugate  transform  of  an  extreme  value  variate,  ^ 
is  the  gamma  variate  of  parameter  w,  and  N(0,1)  is  the  standard  normal 
variate. 

The  explicit  numerical  evaluation  of  the  distribution  of  a  product 
of  independent  betas  arising  in  multivariate  analysis  under  the  normality 
assumption  [2],  [12]  can  be  expedited  with  the  help  of  the  corresponding 
variate  Uj  =  -log  ,  whose  sums  map  into  the  desired  product.  The  inde¬ 
pendent  sum  needed  is  a  multifold  convolution  which  can  be  performed  with 
speed  and  accuracy  by  the  Laguerre  transform  method  [7],  [8]. 


„  s 


§1 .  A  basic  decomposition  of  beta  variates  and  associated  infinite 
factorability 

The  principal  objective  of  this  section  is  the  following  theorem. 
Theorem  1.1 

JJ  =  -log  £  is  infinitely  divisible  for  any  j3  variate.  Equivalently, 

any  8  is  infinitely  factorable,  i.e.,  6  =  6  *6  _...6  where  <5  .  are 

K  *  ~nl  ~n2  ~nn  ~nj 

i. i.d. 

We  will  prove  this  theorem  through  a  lemma  which  provides  insight  into 
the  structure  of  beta  variates  and  is  of  some  interest  in  its  own 
right. 

Lemma  1 . 2 

Let  w  =  [w]  +  0,  where  [a]  is  the  largest  integer  less  than  or  equal 
to  a  and  0  s  0  s  1.  Then  one  has  the  decomposition 


„  ..CM  ..PF 

(1.1)  U  =  U  -  +  U  „  r  > 

~v,w  ~v,0  ~v+0,[wj 

where  (a)  and  [w]  are  independent;  (b)  U^'*0  has  a  completely 

monotone  p.d.f.  when  0  <  0  <  1  and  =  0.  Furthermore,  U^1*  +  -  E 

r  ~v,0  ~  ~v,0  v  ~ 

PF 

as  0  -*•  1,  where  E  is  an  exponential  variate  with  mean  one;  (c)  Q  ^ 

PF 

has  a  PFw  p.d.f.  when  [w]  =  1,2,3...  and  Uv+0  Q  =  0. 


Proof 


From  (0.3),  one  has 


„  _  fr(v»e)  r(v+s)  Hr(v+e+fwi)  r(v+e+s)  , 

ii.2)  -  i  r(v)  r(v+0+S)n  r(v+0)  r(v+e+fwl+sj; 


=  *li;v,e(s)’^!i;v+0,[w] 


(s)  , 


and 


i.e.,  U 


Vw  s  ^v,e  +  SJve.M  and  V*  =  £vfe%+e,[w]*  where  Ve 


+  u 


and 


liv+ejw]  are  independent  and  j^Q  and  £v+0  ^  are  independent.  The 
density  of  =  -log  g^0  is,  from  (0.2),  fU;v>0(y)  =  e'vy(l  -  e"y)6_1. 

Consequently, 


(1.3) 


1 


fU;v,0(y)  =  B(v,0)  JQ  p0ke 


- (k+v)y 


where  p 


0O  =  1  and  p0k  =  H  [1  -  (0/j)],  0  <  0  <  1,  k  >  1,  so  that 


3  =  1 


.CM 


0(y)  is  completely  monotone.  We  write  U  .  =  IT1*  .  We  note  from 
jd;v,0  r  7  ~v,0  ~v,0 


(1.2)  that,  for  Re  s  >  -v,  (s)  +  1  as  0  +  0  and  therefore  =  0. 

u,v,  e  ~v,  o  ~ 


Similarly,  Q(s)  -►  as  0  -*•  1  and  U^M0  +  as  0  +  1,  proving  (a) 


U;v,0 


s+v 


M-i 


and  (b).  For  (c) ,  we  see  that,  for  [w]  £  1,  <j>u>v+0  ^  (s)  =  n  (v+0+j)/ 


(s+v+0+j).  For  the  variate  Uv+0  =  -log  ,  one  has  therefore 


(1.4) 


M-l 

u  =  z 

~v+6.[w]  j _0  v+0+j  ~3 


1 


-  E. 


where  the  Ej  are  independent  exponential  variates  with  E[E^]  =  1.  It 

follows  that  j-wj  has  PF  density  [4],  when  [w]  >  1,  and  we  write 

PF  PF 

Uy+0>[wj.  From  (1.2)  we  see  that  JJy+Q  Q  =  0,  proving  the  lemma.  □ 

Proof  of  Theorem  1 . 1 

It  has  been  shown  by  F.  Steutel  [13]  that  any  completely  monotone 
variate  is  infinitely  divisible.  Since  E  is  infinitely  divisible  and 
the  sum  of  infinitely  divisible  variates  is  infinitely  divisible,  the 
result  is  immediate.  □ 


The  decomposition  of  Lemma  1.2  shows  that  f^(x)  is  the  convolu- 

/*w 

tion  of  a  strongly  unimodal  p.d.f.  [3],  the  PF  component,  and  a  com¬ 
pletely  monotone  component,  shedding  additional  light  on  the  familiar 
unimodality  of  all  beta  variates. 

Remark  1 . 3 

As  shown  in  [9],  any  p.d.f.  f^(x)  with  the  decomposition  (1.1) 
has  the  property  that  f^(x)*f  ^(x),  where  the  asterisk  denotes  convo- 
lution,  is  a  scale  mixture  of  symmetric  normals,  and  that  for  such  a 
distribution,  distance  to  normality  is  measured  by  the  kurtosis  of  X. 

The  kurtosis  of  U,  therefore,  provides  a  consistent  measure  of  the  log- 
normality  of  £  described  in  the  next  section. 

The  decomposition  (1.1)  has  also  been  demonstrated  in  [6]  for  any 
passage  time  Tmn  between  any  two  states  m,  n  of  any  birth-death  process. 


■  « 


s  * 


i 


. 
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§ 2.  Asymptotic  behavior  of  beta  variates  for  large  v  and  w 

We  turn  next  to  the  asymptotic  behavior  of  the  U  and  Jj5  variates. 
Theorem  2. 1 

Let  w  =  [w]  +  0,  0  £  0  <  1.  Then  for  any  v  >  0,  Ur  -  log  w  ^  G 
as  w  +  +«  where  the  p.d.f.  of  G  is  given  by  f_(v)  =  (e  *-v  *^/T(v))* 

^  o' 

-y  ~ 

-y  -e 

e  -oo  <  y  <  oo 


Proof 

PF 

Let  K  =  [w]  -  1  so  that  w  =  K+1+0.  Let  S  =  U  .  r  ,  -  log  w. 

1  J  ~v,w  ~v+0,[w] 

Then,  from  (1.4),  the  Laplace  transform  of  the  p.d.f.  of  S  is  given 
K  ~V’" 

by  <<>g.v  (s)  =  n  (v+0+j)  (K+1+0) S/ (s+v+0+j )  .  This  can  be  rewritten 

’  j=0 

as 


v+0 


i  ..s+v+0, 
~rrrr  ^  t  • 


(2.1)  ^S;v,w^s^  ~  s+v+B  s+v+0+j 


{  n  k'(v+0)}-(i  +l^)s  . 

j=i  3  K 


The  first  bracket  in  (2.1)  converges  to  T(l+s+v+0)  while  the  second 
converges  to  l/T(l+v+0)  as  K  +°°.  Hence,  for  0  <  0  <  1  fixed,  one  has 


(2.2) 


*S;v,w(s) 


r  (s+v+9) 

r(v*6) 


as  w  +<* 


_  ,,  ,,  ,  ,  ,,  T(v+0)  T(v+s) 

From  (1.2)  one  has  $,j.v>e(s)  =  T^T  ‘  f~v^9^s~)  ‘  Slnce 

U  -  log  w  =  U™  +  UPF„  r  T  -  log  w  =  UCMa  +  S  ,  one  has,  from 
~v , w  K  »v,0  *v+0,[w]  5  ~v,0  ~v,w’ 

(2.2)  for  every  fixed  6  (0  £  0  <  1), 


*U;v,w(s:)  =  <t’U;v,0(s-)*'*’S;v,w(s)  ^  T(v) 


U>v,0  "'Sjv.w' 


(2.3) 


,  as  w  +  +•  , 


where  this  limit  is  independent  of  6.  □ 


We  note  that  the  extreme  value  distribution,  whose  p.d.f.  is 
-y  -e'y 

fL(y)  =  e  ,e  ,  -°°  <  y  <  °°,  has  the  bilateral  Laplace  transform 

(s)  =  r(l+s).  Hence  the  limit  of  both  <j>c  (s)  in  (2.2)  and 
A;  b;v,w 

(s)  in  (2.3)  are  conjugate  transforms  [5]  of  <*>.(s). 

v ,  w  L 

Theorem  2.1  above  describes  the  asymptotic  convergence  in  dis¬ 
tribution  of  U  -  log  w  for  v  fixed  as  w  -*•  +<*>,  to  an  extreme  value 
variate.  In  the  next  theorem,  we  deal  with  the  asymptotic  behavior 

for  w  fixed  as  v  +  +«,  and  show  convergence  in  distribution  of  vll 

'"W,V  ,  w 

to  a  Gamma  variate.  Finally,  in  Theorem  2.5  we  will  be  dealing  with 

sequences  (v  ,  w  )  in  which  both  v  and  w  become  infinite  in  a  spe- 
n  n'  n  n  r 

cified  way,  and  asymptotic  normality  will  be  demonstrated.  The  three 
cases  are  shown  graphically  in  Fig.  1  (a),  (b) ,  (c) . 

Theorem  2.2 

Let  w  =  [w]  +  0  >  0  be  fixed  where  0  <  0  <  1.  Then  one  has 
vUy  w  ^  Xy,»  as  v  -*•  +«,  where  ^  is  the  gamma  variate  with  p.d.f. 
fv  w(s)  =  xW_1e'X/r(w). 

Proof 


L«  C  be  the  P-d-f-  °f 

fCM  ,  }  =  — j — 

U;v,0m  B(v, 0) 

given  by 


e"v^(l  -  e-^-'1-65 


Then,  from  (0.2), 
The  p.d.f.  of  vUCfl' 

r  ^\r  H 


is  then 


I  fCM  (L\  -  I  r(v+8)  .  e  7 

v  U;v,0V  "  v  r(v)T(0)  *  n  .  -y/vJ-9 


-8- 


From  the  Stirling  formula,  T(v)  ~  /2tT  vV'^e  V  as  v  -*■  +»,  one  has 


IfcM  (iwJLye-i 
v  u;v,eV  r (6)  y 


e  '  as  v  -►  , 


CM  d  P  F 

i.e.,  vll  +  y.  as  v  ■>  ®,  For  the  PF  part  U  „  r  ,  with  [wl  >  1,  one 
~v,9  r  ~v+8,[w]  1  J 

has  from  (1.4)  that 


PF  M'1  1  d 

vU  a  r  ,  =  I  - V-r  E.  °  yr  i  . 

~v+0>M  i=1  ,  +  lii  XM 

J  V 

The  theorem  now  follows  from  Lemma  1.2.  □ 

From  (0.3),  one  has 


'2-4a>  V.w  *  ■ 


a? log  V,v,»(s)  s=o '  - <,(v)  • 


(2-4b)  °U;v,w  ■  Var!!iv,J  ■  (3T)2  108  *U;v,wts)ls=.0  *  *’  (v)  *  (v*“)  • 


where 


(2.S)  *(z)  =  ^  log  T(z)  =  r*(z)/r(z)  ,  Re  z  >  0 


(2.6)  Z  =  (U  -  y„  )/o„ 
v  '  ~v,w  V~v,w  U;v,w  U;v,w 


We  next  show  that  Z  -*■  N(0,1)  as  v  and  w  go  to  +»  along  certain  paths. 
~v ,  w 

Two  preliminary  lemmas  are  needed. 

Lemma  2.3 


<  1+x  for  all  x  >  0  . 


1  -  e 


-9- 


Proof 


<*  k 
r  x  x 


It  is  clear  that  1+x  <  £  ry  =  ex  for  all  x  >  0.  Then  (l+x)e  <  1 

k=0 

so  that  x  <  (1+x) (1  -  e~x) ,  and  the  result  follows.  □ 

Lemma  2 . 4 

Let  P  be  any  directed  path  in  the  (v,w)  plane  for  which  either 

(a)  v,w  •+  +®  and  v  <  Kwa,  0  <  a  <  1,  K  >  0,  or  (b)  v,w  -+  +®  and 

v  2  Kw  ,  a  2  1,  K>0.  Then  va„  ■+  +»  as  v,w  ■+  +°°  along  P. 

u;v  ,w 

Proof 

We  note  the  following  identity  [1] 


i  "  t  _zt 

(2.7)  ip  (z)  =  /  — - r  dt  ,  Re  z  >  0  . 

0  1  -  e 


OQ 

Then,  from  (2.4b),  v2a.2  =  v2  j  - - — r  e"vt(l  - 

-t  ’  0  1  -  e"1 

since  - - —  <  1  for  t  >  0,  we  have 


-Wf 

e  )dt.  Hence, 


w 

00  -  -  X 

v2a2  2  v  /  e'X(l  -  e  V  )dx 

U;v>w  J 


a  1  -a 

Ifv<Kw,0<a<l,  K  >  0,  one  has  w/v  2  w  /K  and  by  the  dominated 

•  “X  -wx/v 

convergence  theorem  /  e  (1  -  e  )dx  -+  1  along  such  a  path.  For  the 

0  1-a 

case  (b) ,  one  sees  that  w/v  sw  /K,  a  M,  K  >  0  and  w/v  ■+  0  along  P. 
One  then  has 


v  o.,  2  K  w  /  e 

u;v»w  0 


■x  1  -  e 


w 

-  —  X 
V 


dx 


The  integral  in  the  last  term  converges  to  /  xe"  dx  =  1  by  the  dominated 

0 

convergence  theorem  and  the  lemma  follows.  □ 


-10- 


We  are  now  ready  to  show  asymptotic  normality  of  U  under  the  condi- 

,  w 

tions  of  Lemma  2.4. 

Theorem  2.5 

Z  -*■  N (0 , 1 )  as  v  and  w  go  to  +»  along  any  path  P  as  given  in 
~v ^ 

Lemma  2.4. 

Proof 


We  write  p  =  p..  and  o  =  a..  for  notational  simplicity.  It 
U;v,w  U;v,w  1 

is  clear  from  (2.6)  that 


sZ  — 

C2.8)  *  (s)  =  E[e  ~]  -  e7  \.  (|) 

'  (V'  ' 


We  note  that  for  sufficiently  small 
v+w  v+w+s 


Vv  wW  =  exp[  /  *(u)du  - 

»  v 

v+w+s 

/  tj;(u)du]  so  that 


/ 

v+s 


|s|,  one 
i|)(u)du]  * 


v+w 


has  from  (0.3)  and 
v+s 

exp[  /  ij>  (u)du  - 
v 


(2.5) 


(2.9) 


s 

4’,,..,  =  exp[/  {'Kv+X)  -  'l'  (v+w+x)  }dx] 

x»v»w  0 


<J>_.  (s)  in  (2.8)  can  then  be  rewritten  from  (2.9)  as  4>_  (s)  = 

j  w 

s/a 

exp[  /  {(iKv+x)  -  *(v))  -  (i|» (v+w+x)  -  il<(v+w))}dx]  .  By  letting 

0 

y  =  ax,  we  obtain 

s 

/  h(v,w,y)dy 
(2.K0  f  (s)  =  e° 

<>s» 7  ' 

where 


1  Y,r°  1 

h(v,w,y)  =  —  /  (v+u)  -  (v+w+u))du 

®  A 


(2.11) 


-li¬ 


lt  will  be  seen  that  ^  h(v,w,y)  -*•  1  for  all  y  >  0  as  v,w  -►  +®  along 

1  ^ 

the  path  given  and  that  one  then  has  <J>  (s)  -*■  e^5  ,  as  needed. 

>  v  >  w 

From  (2.11),  h(v,w,y)  =  (v  +  y/a)  -  (v  +  w  +  y/o) }  so 

that  (2.4b)  and  (2.7)  lead  to 


(2.12)  ^  h(v,w,y)  -  S-U-S 


/  — *-(v  4  y/<,)t(l  -  e~Kt)dt 

n  1  *  “  t 


r  t  ,-vt„  -wt, 
i  - r  e  (1  -  e  )dt 

0  1  -  e 


We  note  that  ^  h(v,w,y)  is  monotone  decreasing  in  y  (y  >  0),  and 

0  S  5y  h(v’w«>r)  ~  1  for  a11  v>w>y  >  0.  Let  x  =  vt.  Then  (2.12) 
becomes 


(2.13)  hCv,w,y)  = 

From  Lemma  2.3,  - *^V- 7—  <  1  +  ^  for  x,v  >  0  and  by  the  dominated  con- 

1  .  e”x/v  v 

vergence  theorem  one  can  pass  v  to  the  limit  along  the  path  given.  It 
follows  from  Lemma  2.4  that  h(v,w,y)  +  1  as  v,w  -+  +®  along  the  path 

y.  d 

given.  From  (2.11),  h(v,w,0)  =  0  so  that  h(v,w,y)  =  /  ^  h(v,w,u)du. 

d  0 

Since  0  s  ^  h(v,w,y)  <,  1,  one  sees  that  h(v,w,y)  -*•  y  as  v,w  +®  along 

the  path,  again  by  the  denominated  convergence  theorem,  for  any  y  >  0. 

The  theorem  then  follows.  0 

The  convergence  vU  described  in  Theorem  2 . 2  has  been  shown  by 
G.  S.  Mudholkar  and  M.  C.  Trivedi  (private  communication).  They  also 


/  - — 

0  1  -  e 


x/v 

-x/v 


-O  +  y/ov)x  -  ~  x 

e  (1  -  e  )dx 


0  1  -  e 


x/v 

-x/v 


-X,, 

e  (1  -  e 
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state  that  U  is  "asymptotically  normal  as  v,w  +»"  but  do  not 
~v ,  w 

provide  a  proof  [10]. 

In  the  original  form  of  Theorem  2.5,  only  ray  paths  v  =  Kw,  K  >  0, 
were  considered.  A  referee  suggested  the  more  general  paths  of  Lemma 
2.4,  and  indicated  that  the  result  might  also  be  obtained  from  Chap.  4, 
Theorem  18  of  V.  V.  Petrov  [11]. 


U  -  log  w  +  G  as  k  +  ♦»  along  P. 

~V  y  W  ^  A 


vU  ^  y  as  v  -*•  +®  along  P- 

'"'V  y  W 


z 

~v,w 


u 

~v,w 


V  ,  ' 


0,, 

U;v,w 


N(0,1) 


as  v,w  -*■  +<»>  along  P^,  where  P^  is  a  path  such  that  either 
0  <  a  <  1 ,  K  >  0  or  v  i  Kwa ,  a  2  1 ,  K  >  0 . 


Fig.  1.  Asymptotic  Behavior  of  Uy  w  =  -log  £y  w 
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§3.  Generation  of 


random  numbers 


CM  PF 

The  decomposition  U  =  U  +  u  r  ,  in  Lemma  1.2  may  be 

f  W  |  o  '"'V  ^  0  j  I  W  I 

employed  to  provide  a  simple  algorithm  for  generating  8  .  random  num- 

,  w 

CM 

bers.  From  (1.3)  the  Laplace  transform  of  the  p.d.f.  of  ^  can  be 
given  by 


<3-15  W.9(s)  ■  i  «k 

v  =  n 


where 


(3.2) 


’k  "  B(v, 0) (v+k)  *  6k 


_  (eC1)(-l)k  ,  k  =  0,1,2,...  . 


*  0  K 

It  is  clear  that  q,  >  0  for  all  k.  One  sees  quickly  that  £  — r  = 


k=0  *  v  *  0 


•  pev 

(3.3)  l  ~--B(v,0) 

k=0  v  1 


and  therefore  (q^)g  is  a  probability  distribution.  Let  EL  be  i.i.d. 
with  the  common  c.d.f.  1  -  e  X,  j  =  0,1, — ,M  =  [w] .  From  Lemma  1.2, 
(1.4)  and  (3.1),  one  then  has 


(3'4)  =  j0  ^  h 


where  N  is  the  discrete  random  variable  with  P[N  =  k]  *  and  independent 
of  E...  Let  U.  be  independent  and  identical  uniform  variates  on  (0,1).  Since 
Uj  ■  e  i  and  U,>w  a  -i°g  &vtW>  Etl-  (3.4)  leads  to 
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(3.5) 


&v,w 


M-l 

n 

3=0 


1 

yV+0+ j 

~j 


1_ 

N 


Hence  one  has  the  following  algorithm  for  generating  w  random  numbers. 
Algorithm 

(a)  Generate  [w]  +  l  independent  and  identical  uniform  variates  LL(u)), 
j  =  0,1,... ,M  =  [w],  on  (0,1). 

oo 

(b)  Generate  the  variate  N(uj)  from  the  distribution  (qk)Q- 

M-l  — N(w) 

(O  w(“)  =  n  yi(“)V  J  ' 

^  ’  j  =0  J 

The  algorithm  is  simple  and  straightforward.  Advantages  and  disadvantages 
of  the  algorithm  with  respect  to  existing  algorithms  will  be  described 


elsewhere. 
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§4.  Explicit  calculation  of  the  distribution  of  the  product  of  independent 
Beta  variates  ™ 

For  certain  likelihood  ratio  statistics  arising  in  multivariate 

analysis,  one  must  evaluate  the  distribution  of 


(4.1) 


X  = 


£v 


w  %  w  ‘ 
1*1  2*  2 


K,WK 


where  the  beta  variates  are  independent.  This  distribution  may  be 
obtained  via  the  Laguerre  transform  procedure  described  in  [7],  [8]  in  the 
following  way.  From  (4.1) 

K  K 

(4.2)  -log  X  =  l  (-log  w  )  «  l  Uy  w  • 

3  =  1  3*  3  3  =  1  y  3 

The  U  variates  are  independent  and  absolutely  continuous  with 

~Vj,W. 

p.d.f.'s  as  in  (0.2).  They  therefore  have  the  properties  of  regularity 
and  rapid  decrease  required  by  the  Laguerre  transform  method  for  con¬ 
volving  p.d.f.'s  and  permit  vector  representations  of  modest  length 
with  high  accuracy.  The  Laguerre  transform  coefficients  required  are 
easily  obtained  analytically  and  the  calculation  of  the  p.d.f.  of 
-log  X  and  hence  of  X  proceeds  rapidly. 
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